Sound emission due to superfluid vortex reconnections 
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By performing numerical simulations of superfluid vortex 
ring collisions we make direct quantitative measurements of 
the sound energy released due to vortex reconnections. We 
show that the energy radiated expressed in terms of the loss 
of vortex line length is a simple function of the reconnection 
angle. In addition, we study the temporal and spatial distri- 
bution of the radiation and show that energy is emitted in the 
form of a rarefaction pulse. The pulse evolves into a sound 
wave with a wavelength of 6 — 8 healing lengths. 
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The decay of superfluid turbulence in the limit of low 
temperature raises some fundamental questions in the 
field of quantum fluids. For example, vortex tangles pro- 
duced in superfluid helium at T < 100 niK are observed 
to decay m , but at these temperatures the frictional dis- 
sipation due to thermal excitations is practically negli- 
gible. Recent theoretical work [0-0 has highlighted the 
possible role of sound emission as a dissipation mecha- 
nism. Sound emission may occur due to vortex motion 
or reconnection, however, conventional numerical simula- 
tions based on vortex filaments governed by incompress- 
ible Euler dynamics (the Biot-Savart law) are unable to 
describe this process. A useful tool to study vortex- 
sound interactions is the Gross-Pitaevskii (GP) equation. 
Although unable to fully represent the physics of Hell, 
the GP equation does provide a sophisticated fluid dy- 
namical model capable of describing vortex nucleation 
B and reconnections &-@l. Furthermore, the GP equa- 
tion has been shown to provide an accurate description of 
the recently discovered atomic Bose-Einstein condensates 
PI , where similar issues of vortex-sound interactions are 
likely to occur, particularly in experiments where many 
vortices are formed [lOUO]. 

In this letter, we compute the sound energy radiated 
due to superfluid vortex reconnections using the GP 
model. We parameterise the energy in terms of vortex 
line length, and show that the vortex line length de- 
stroyed during a reconnection is a simple function of the 
reconnection angle. In addition, we measure the tem- 
poral and spatial distribution of the radiation and show 
that energy is emitted in the form of a rarefaction pulse 
which subsequently disperses into sound waves. 

To produce a reconnection we collide two vortex rings 
whose axes of propagation are offset. This system has the 
advantage that by varying the offset we can study a range 



of reconnection angles. In addition, a well behaved ini- 
tial state can be constructed within a spatially conflned 
region. We adopt the same numerical methods as in our 
previous work [p^-p^ . The key points are that through- 
out the paper we use dimensionless units, where distance 
and velocity are measured in terms of the healing length, 
^, and the sound speed, c, respectively. In addition, the 
asymptotic number density, no, is rescaled to unity. The 
initial state is taken as the product of two vortex ring 
states, il}{r,Q) — (j)i{r — ri)(j)2{r — r2), where 0i,2('") 
are time- independent vortex ring solutions of the uniform 
flow equation found by Newton's method |13|. We ar- 
range for the rings to propagate in the ±x directions with 
an offset in the y direction, i.e., ri,2 = (±25u, ibD/2,0), 
where a;, y, and z correspond to the axes shown in Fig. 1, 
V is the ring velocity, and D is the offset between the 
propagation axes of the two rings. The initial state is 
evolved according to the dimensionless GP equation. 
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using a semi-implicit Crank-Nicholson time step. To 
model a large box, we map an infinite length onto the 
space — 1 < a;' < 1 using x' — x/{\x\ + Q with C = 12. 
We use a grid spacing. Ax' = 0.007, and a time-step. 
At = 0.01. Simulations have been performed with vortex 
ring radii of i? = 5.04, 6.00 and 7.30, which correspond 
to velocties v — 0.34, 0.3, and 0.26, respectively. 

A typical sequence illustrating the vortex ring colli- 
sion is shown in Fig. 1. As the rings approach they 
stretch in the yz plane. The reconnections occur along 
the z axis at around i = 30. The reconnections produce 
two highly elongated rings and two sound pulses which 
propagate outwards along the ±z axis. The sound pulse 
appears in Fig. 1 as an oval in the centre of the top view 
at t = 40. The stretched rings (t = 50) rapidly shrink 
into two smaller vibrating rings which move outwards. 
For the offset 13 = 4 shown in Fig. 1, the outgoing rings 
propagate at an angle (j) — ±56° to the x axis. The 
smaller radius of the outgoing rings reflects the energy 
loss due to the reconnections. For larger offsets, the col- 
lision is less violent than the example shown, and the 
scattering angle of the outgoing rings, 0, is smaller. 

To quantify the sound energy we calculate the energy 
within a measurement sphere of radius 20. The energy 
is defined relative to a uniform laminar state as in our 
previous work p3| . The energy per ring as a function 
of time for collisions between rings with radii R — 6.0, 
offset hy D — 2 to D — 8.5, is shown in Fig. 2(inset). If 



the reconnection occurs a few healing lengths along the 
z axis at i '--^ 30, and the sound pulse travels outwards at 
speed c—1, then one would expect to observe the sound 
leaving the measurement volume between i = 45 and 
t = 50. This appears as an energy drop in the numerical 
data, Fig. 2(inset). After the sound pulse has left the 
sphere one expects the energy to remain approximately 
constant until the outgoing vortex rings leave between 
t = 60 and i = 75. This energy 'plateau' is well defined 
in Fig. 2 (inset) for offsets between D = A and D = 8.5, 
but for smaller offsets the outgoing rings are too fast to 
be resolved from the sound pulse. We define the radiated 
sound energy as the difference between the initial energy 
and the value at the centre of the plateau. 
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FIG. 1. Sequence of density isosurfaces {\ip\^ = 0.75) illus- 
trating a vortex ring collision with ring radius R — 6 and 
offset D = 4. A side view (a) and a top view (b) are shown. 
Initially, the rings propagate along the x axis and collide in the 
yz plane at £ ~ 30. The collision produces highly stretched 
rings (t = 50) which snap back into two smaller rings moving 
at ±56° to the x axis. The sound pulse is emitted along the 
±z axis, and appears as an oval in the centre of the top view 
at t = 40. 

A convenient way to express the energy loss is in terms 
of an effective vortex line length destroyed. To calculate 
the change in vortex line length we convert between en- 
ergy and length using the analytical result for the energy 
of a vortex ring with radius R |Tq] , 



E = 27t^R 
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with a = l/\/2- Clearly, this is approximate as the out- 
going vortex rings are excited. The vortex line length 
destroyed, hoss, as a function of reconnection angle, 
9 = 2cos^^{D/2R), is plotted in Fig. 2. For intermediate 
angles, the data appear to lie on a unique curve, inde- 
pendent of the initial ring radius. A good fit is given by, 
hoss = tan2(6'/2) = (4i?2 - D'^)/D^. For large 9 (smah 
offset D) most of the energy is emitted as sound and 
the energy plateau is not well defined leading to a large 
error. Additional errors occur when comparing Eq. (0) 
with our numerical energy data due to the finite size of 
our measurement volume. This error is largest when the 
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incoming and outgoing rings are largest, i.e. for large R 
and large D. This leads to the over estimate of vortex 
line loss apparent for the R — 7.30, 9 < 7r/2 data in 
Fig. 2. 
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FIG. 2. Vortex line length destroyed as a function of the 
reconnection angle. The data correspond to ring radii of 
R = 5.04 (o), 6.00 (D), and 7.30 (x). The vortex line 
length destroyed is calculated from the energy decrease within 
a measurement sphere when the sound pulse exits. The 
time-dependence of the energy per ring within the measure- 
ment sphere for ring radii R = 6.0 and offsets D = 2 — 8.5 
in increments of 0.5 is shown inset. The sound energy exits 
between t = 48 and i = 52. The energy plateau at i '^ 60 
indicates the energy loss due to sound emission. 

The dependence of the vortex line loss on reconnec- 
tion angle can in principle be applied to determine the 
sound energy emitted in other reconnection geometries. 
However, in the case of straight line vortices, it is known 
that the lines tend to become more anti-parallel as they 
approach |l^,|l^. This will tend to increase the energy 
radiated. To convert the dimensionless units into val- 
ues applicable to Hell, we take the number density as 
TT-o — 2.18 X 10^* m~^, the quantum of circulation as 
K = h/ra = 9.98 x 10^^ m^s^^, and the healing length 
as S,/V2 = 0.128 nm |1^, then the unit of energy is 
Tiuqc^'^ = 6.56 X 10"^** J or 0.475 K. For these parame- 
ters, the energy radiated for reconnection angles of 7r/2 
and 37r/4 are 9.2 x lO'^s J (6.65 K) and 5.2 x lO^^z j 
(37.5 K), respectively. 

To determine the character of the sound radiation we 
have studied the temporal and spatial distribution of the 
emission. The density along the z-axis for a collision 
with ring radius R = 6 and offset 13 = 8 is shown in 
Fig. 3. Initially the density is uniform except for a slight 
increase near the origin indicating the approaching rings. 
Between i = 25 and f = 30 the rings collide in the xy 



plane. During the collision, the vortex cores merge while 
the rings grow outwards. Due to this stretching effect, 
the reconnection point is pushed outwards along the z 
axis to z ~ ±7 and delayed until t ^ 29 (if there were 
no distortion of the incoming rings, one would expect the 
reconnection to occur at z — ±4.5 and t — 25). When 
the circulation is cancelled at i = 29, the density is zero 
at z = ±7. This density minimum continues to move 
outwards in the same direction as the vortex cores prior 
to the reconnection. As the pulse moves outwards the 
density mimimum gradually fills in. 




FIG. 3. The density along the z axis for a collision between 
two vortex rings with radius R = 6 and offset D = 8. The 
eleven curves, corresponding to times i = to i = 50 in incre- 
ments of 5, are offset in increments of 0.25 along the y axis. 
Initially the density is uniform except for a slight increase near 
the origin indicating the approaching rings. Between t — 25 
and t = 30 the rings collide in the xy plane. This results in 
the formation of a rarefaction pulse which moves outwards. 
Note that initially the density at the center of the rarefaction 
pulse is zero and then increases. For t > 40 the pulse evolves 
into a sound wave with a wavelength of 6 — 8 healing lengths. 

As shown for the case of vortex nucleation by a mov- 
ing sphere |]l4| , further insight can be gained by consid- 
ering the time evolution as a transition between time- 
independent states. For the head-on collision {D = 0), 
one can regard the segments of the ring as anti-parallel 
vortex lines. The two dimensional time-independent so- 
lutions for two opposite sign vortices have been studied 
by Jones and Roberts 1 16|. They show that when the vor- 
tex cores merge they form a rarefaction pulse. In a three 
dimensional situation, as the rarefaction pulse expands 
outwards the energy per unit length decreases which cor- 
responds to a lower energy on the dispersion curve. The 
lower energy rarefaction pulse has a higher central den- 
sity and moves faster eventually approaching the sound 
speed |lq| . This exactly describes the behaviour apparent 



in Fig. 3. Eventually the rarefaction pulse evolves into a 
sound wave with a central wavelength of approximately 
7 heahng lengths. For Hell, taking £,/y/2 = 0.128 nm 
||l9| this converts to 1.3 nm, which corresponds to an 
intermediate phonon wavelength (cf. maxon and roton 
wavelengths of 0.6 nm and 0.3 nm, respectively) 

The spatial extent of the sound pulse in the x = 
plane for the same parameters as in Fig. 3 but at a later 
time, t = 54, is shown in Fig. 4. The density variation 
around the origin is a remnant of the departing vortex 
rings. The sound pulse appears as two crescent shaped 
density waves with an angular distribution roughly equal 
to the reconnection angle (indicated by the white lines in 
Fig. 4). A similar distribution is observed in the j/ = 
plane. By t = 54, the sound energy is spread over a large 
area and the amplitude of the density variation is very 
small (only a few percent). The spherical wave fronts 
are consistent with the reconnection position and time 
indicated by Fig. 3. 
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FIG. 4. Density cross-sections in the a; = plane at f = 54 
for ring radii R = 6 with offset D = 8 (same parameters as in 
Fig. 3). The sound pulse appears as two arcs with radius of 
curvature, 25, suggesting that the reconnections occurred at 
z — ±7 and i = 29 consistent with Fig. 3. The angular spread 
is approximately equal to the reconnection angle {6 = 96° for 
this example). The white lines indicate the positions of the 
reconnections, z = ±7, the reconnection angle, 9, and the 
expected position of the sound pulse. Grey scale: black 0.95; 
white 1.025. 

In summary, we have made direct quantitative mea- 
surements of the sound energy emitted due to superfluid 
vortex reconnections. We show that the energy radiated 
can be parameterised in terms of the loss of vortex line 
length which is a simple function of the reconnection an- 
gle. The energy radiated increases dramatically as the 
lines become more anti-parallel which suggests that re- 



connections may be a significant decay mechanism for su- 
perfluid turbulence in the limit of low temperature. We 
also show that the radiation initially appears in the form 
of a rarefaction pulse which evolves into a sound wave 
with a wavelength of 6 — 8 healing lengths. 
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